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LAURA DE MARCO AND AARON SCHIFF 

Abstract. We provide an algorithm for computing the Euler characteristic of the 
curves Sp in "PI™ ~ , consisting of all polynomials with a periodic critical point of 
period p in the space of critically-marked, complex, cubic polynomials. The curves 
were introduced in |Mii IBKMj , and the algorithm applies the main results of }DPj . 
The output is shown for periods p < 26. 



1. Introduction 

Let V^"^ denote the space of cubic polynomials with marked critical points. It 
is convenient to parametrize the space V^"^ by {a,v) G C^, where the pair {a,v) 
corresponds to the polynomial 

fa,v{^) = - 'ia^z + 2a^ + 

with critical points at ±a and critical value v = fa,v{+0')- 

In this article, we study the geometry of the curves Sp C Pf™", introduced by 
J. Milnor in [Mi], consisting of cubic polynomials fa,v for which the critical point +a 
has period exactly p. That is, 

Sp = {(a, v) G rr : flM = a, («) ^ a for all 1 < A; < p}. 

The curve Sp is smooth for all p [ML Theorem 5.1]. As a (possibly disconnected) 
Riemann surface, the curve Sp has finite type: it is obtained from a compact Riemann 
surface Sp by removing finitely many points. The punctures lie at infinity in the space 
-pcm_ date, the irreducibility of Sp is unknown, though it is shown in |BKMt §8] 
for periods p < 4. 

The goal of this article is to explain an algorithm to compute the Euler character- 
istic of the compactification Sp. In |BKMt Theorem 7.2], it is shown to satisfy 

(1.1) x(Sp) = dp{2-p) + Np. 

The number dp is the degree of the curve Sp, and it is easily computable from the 
defining equation. The number Np denotes the number of ends of Sp, the punctures 
Sp \ Sp. Our contribution is the algorithmic process to compute Np, applying the 
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main results of |DP] . The Euler characteristic of Sp is shown in Table [H to period 
p = 26. 

We remark that the computation of the Euler characteristic xi^p) cannot be han- 
dled by traditional methods beyond the small periods. A quick genus computation 
with Maple™, for example, yielded Euler characteristics for p < 4 and failed to pro- 
vide an output for p = 5 where 1S5 is a curve of degree 80. The degree of Sp is 
on the order of 3^~^, and the curves Sp will be highly singular at infinity for any 
choice of projective compactification of V^"^ ^ and p sufficiently large. The Euler 
characteristics for periods p < 4 appear in |BKM] . 

1.1. Outline of the algorithm. As described in |Mi] . the ends of Sp correspond 
to the escape regions of Sp, the open subsets of Sp consisting of polynomials with 
the critical point —a tending to infinity under iteration. The main ingredient in 
the computation of Np is the combinatorial analysis of polynomial dynamics on the 
basin of infinity, developed in |BH] and |DP] . Recall that the basin of infinity of a 
polynomial / is the domain 

Xif) = {zeC: riz) ^ 00}. 

From |BH] . we use the properties of the tableau (or equivalently, the Yoccoz tau- 
function) of a cubic polynomial; this combinatorial object encodes the first-return of 
a critical point to its "critical nest." From [DPj we use the combinatorics of the pic- 
tograph, a more refined encoding of the first-return of a critical point to a "decorated 
critical nest," allowing us to distinguish and count topological conjugacy classes. 
The steps of the algorithm are: 

(1) Fix p. For each k dividing p, with 1 < k < p, determine all admissible 
tau-functions with period k. 

(2) Count the number of topological conjugacy classes of basins of infinity (/, X(/)) 
associated to each tau-function. 

(3) Compute the number of topological conjugacy classes of polynomials in Sp 
with one escaping critical point: each class is determined by the class of its 
basin of infinity (with a tau-function of period k) and a point in the Mandel- 
brot set associated to a period p/k critical point. 

(4) Determine the number Np of escape regions in Sp-. there are either one or 
two ends in Sp associated to each topological conjugacy class computed in the 
previous step, determined by the twist period of the tau-function. 

(5) Test the output against the degree of Sp-. Np is the total number of escape re- 
gions, while the degree of Sp must equal the number of escape regions counted 
with multiplicity. The multiplicity is computed from the tau-function. 

Step (1) uses the tableau rules of |BH] . as corrected in jKi ] IDM] : a translation into the 
language of the Yoccoz tau-functions was given in [DSj . The bulk of the computing 
time and memory usage goes into Step (1). In ^ we provide the theoretical results 
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13639477 
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24 
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25 
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Table 1 . The output of the Euler Characteristic algorithm. From left 
to right: the period p; the number of tau-functions with period p; the 
number of escape regions of Sp with the hybrid class of z'^ (see Theorem 
15. 3p : the Euler characteristic xi^p)'^ and a comparison to 3^~^. 



needed for the computation. We include the theoretical results we used for improving 
the speed of the algorithm; we believe that some of these are interesting in their own 
right. 

Step (2) was implemented already in [DSj . applying the results of |DP] . Step (3) 
relies on the work of Branner and Hubbard in |BH] (see also |BKMt Theorem 3.9]), to 
know that the conformal class of a cubic polynomial in an escape region depends only 
on the class of its basin and the class of its degree 2 polynomial-like restriction. Steps 
(4) and (5) are explained in ^ where we relate an escape region in Sp to its quotient 
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in the moduli space of cubic polynomials M.'^. The multiplicity of an escape region 
is computed and depends only on the underlying tau-function. 

1.2. Details of the computation. An implementation of the algorithm was written 
with C++. We compiled the output in Table [T]to period p = 26. The low periods are 
computed quickly, while the computation for period 26 took 9:13 hours (Intel Core 2 
Quad @ 2.5 GHz on Windows 7 32-bit edition), executed on a single thread. 

1.3. The growth rate of x{<Sp)- An easy computation shows that —x{<Sp) — )■ oo as 
p — i- oo [MiJ. Using methods from pluripotential theory, Dujardin showed that 



as j9 — )■ oo [Du]. After viewing the output of this algorithm, Milnor asked whether 
we have 



We include the ratio -xi<Sp) /3P'^ in Table [H 

1.4. Acknowledgments. We would like to thank Jan Kiwi, Jack Milnor, and Kevin 
Pilgrim for helpful conversations and attention to the output. 



In this section, we define the Yoccoz tau-function of a cubic polynomial and explain 
Step 1 of the algorithm, the procedure to compute all periodic tau-functions of a given 
period p. The main theoretical result is the following: 

Theorem 2.1. For each period p >l, a tau-function has period p if and only if 

T{n) = n — p 

for all n > 2p ~ 2. 

We show that the bound 2p — 2 is optimal: for every p > 3, there exists a (unique) 
period p tau-function with r(2p — 3) p — 3. See Lemma \2M 

As described below, it is quite easy (from a theoretical point of view) to generate the 
periodic tau-functions, combining Theorem 12.11 with Theorem 12.21 A first approach 
might be to generate all admissible tau-functions of length 2p — 2 and test for equality 
t{2p — 2) = p — 2. As witnessed by the computations of [DSj . however, the number 
of tau-functions grows exponentially with length, and only a small proportion are 
periodic. For example, there are 649 tau-functions of period p = 10, while there are 



(1.2) 




as p — )■ oo. O: 



'r, equivalently by equation (11. ip . do we have 

Np = 0(3^"^)? 



2. The r functions 
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279,415 tau-functions of length 2p — 2 — 18. Much of this section is devoted to the 
results we apply to reduce the computation time and memory usage. 

2.1. The tau-function of a polynomial. Fix a cubic polynomial / with discon- 
nected Julia set, and let 

Gf{z)^ hm ^log+|r(^)| 

be its escape rate. Let Ci and C2 be the critical points of /, labeled so that Gf (02) < 
Gf{ci). For each integer n > such that Gf{c2) < G/(ci)/3"~^, we define the 
critical puzzle piece Pn{f) as the connected component of {z : Gf{z) < G'/(ci)/3"'~^} 
containing C2. The puzzle piece Po{f) contains both critical points. For positive 
integers n, we set 

r(n) = max{j < n : /"-■'(cs) G P,(/)}, 

defining a function r from {1, . . . , N} (or all of N) to the non-ncgativc integers. The 
largest on which r is defined is said to be the length of the tau-function. In other 
words, N is the greatest integer such that Gf (02) < G/(ci)/3^~^. If there is no 
maximal N, we say r has length 00. 

The markers of a tau-function with length N are the integers 

{m e {1, . . . , - 1} : r{m + 1) < r{m) + 1}. 

The marked levels of r are all integers in the forward orbit of a marker: 

{I > : I — T'^{m) for marker m and n > 0} U {0}; 

we say is marked even if there are no markers. The positive marked levels coincide 
with the lengths of the columns in the Branner-Hubbard tableau. In terms of the 
polynomial /, a level / is marked if the orbit of the critical point intersects Pi{f) \ 
P;_i(/). We say / is marked by k if the fc-th iterate f^{ci) lies in \ P/_i(/). 

2.2. Properties of tau-functions. Let N denote the positive integers {1, 2, 3, . . .}. 

For any positive integer N , a function 

t:{1,2,3,...,7V}^NU{0} 

or a function 

T : NU{0} 

is said to be admissible if it satisfies the following properties (A)-(E): 

(A) r(l) = 

(B) T(n + 1) < T{n) + 1 

From (A) and (B), it follows that T{n) < n for all n e N; consequently, there exists 
a unique integer ord(n) such that the iterate r°'^'^(") (n) = 0. 

(C) If T(n + 1) < r'^in) + 1 for some < A; < ord(n), then T(n + 1) < T''+\n) + 1. 
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(D) If T{n + 1) < T''{n) + 1 for some < k < ord(?T,), and if T{T^{n) + 1) = 
T^+\n) + 1, then r(n + 1) < T^+\n) + 1. 

(E) If ord(n) > 1 and ord(r°'''^(")-i(n) + 1) = 1, then r(n + 1) ^ 0. 

A tau-function is admissible if and only if it is the tau-function of a cubic poly- 
nomial |DSl Proposition 2.1]. The proof is by induction on N, applying the rules 
for admissible tableaux in [BH] . Property (E) is another formulation of the "missing 
tableau rule" (M4) appearing in [Ki] and [DM]. 

Let k be the number of markers which appear in the orbit 

N ^ t{N) ^ ...^ T^'^^i^^N) = 0, 
and label these k markers by l[, /g, ■ • • , so that 

iv = > /; > > ■ ■ ■ > z;^ > 0. 

For each < i < k, let k = r(Z-) so that 

r{N) = lo>h> ■■■>lk>0. 
Properties (A)-(E) imply the following: 

Theorem 2.2. |DSl Theorem 2.2] Given an admissible tau-function r of length N, 
an extension to length N + 1 is admissible if and only if 

r{N + 1) = /j + 1 for some <i <k 

or t{N + l) = 0iflk>0ork = 0. 

Note, in particular, that t{N + 1) = t{N) + 1 is always an admissible extension to 
length 1. 

2.3. Periodic tau-functions. For cubic polynomials with exactly one critical point 
in the basin of infinity, the tau-function will have infinite length. An admissible tau- 
function r : N — 7- N U {0} is periodic with period p if there exists N{t) G N such 
that 

T{n) = n — p 

for all n > N{t). Such tau-functions correspond to basins of infinity with a bounded 
critical orbit in a periodic component of the filled Julia set; r has period p if and only 
if the component has period exactly p. For computational purposes, we need a bound 
on N{t) depending only on the period p. The bound N{t) < 2p — 2 is granted by 
Theorem 12. ![ which we prove below. 

Lemma 2.3. Ifr has period p, thenT{n) > n — p for alln. Further, ifT{nQ) = n^—p 
for some Uq, then T{n) = n — p for all n > Uq. 

Proof. This follows easily from property (B). □ 
Lemma 2.4. If t has period p, then I < p — 1 for all marked levels I. 
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Proof. Let / be any cubic polynomial with a given periodic tau-function. Label the 
critical points of / as in §2.11 Without loss of generality, we may assume the critical 
point C2 is periodic with period exactly p. 

Suppose I is marked by iterate k, and assume first that r(/) = 0. From Lemma 
I2.3[ we have 1 = 1 — t{1) < p. The first return of Pi to the critical nest occurs with 
f\Pi) = Pq- Because it maps with degree 2, the iterates /'(C2) and f^^\c2) must lie 
in the two distinct components of {Gf < G/(ci)} inside Pq. By periodicity, then, we 
must have I < p. 

More generally, we have that the first return of Pi to the critical nest is f^~'^'^^\Pi) = 
Pt{i), and / — r(/) < p. As above, because the first return is with degree 2, the images 
^'--^(0^2 and f^^^~'^''^\c2) cannot lie in the same component of {Gf < G/(c2)/3^^'^} 
within Pt{i)i while C2 and f^{c2) do lie in the same component. Therefore 1~t{1) < p. 

In addition, we must have / — t^(j9) < p, as this is the first level where the forward 
orbit of C2 and f^{c2) might come together. If / is not a marker, then /'~^*-'^(c2) lies 
in the same component as C2 at t(/), and therefore its image at t^(Z) is in a distinct 
component from that of '■'•'(02). On the other hand, if / is a marker, then r(/) 

is marked by / — r(/). By periodicity, we can take k = p — {I — t{1)). At t^(/), we 
have /'~^^*-'-'(c2) and /^+'~'^^(')(c2) again in distinct components. In either case, we 
conclude that / — r^(/) < p. 

We continue inductively. For the induction step, we begin with I — t'^{p) < p and 
/-r"+i(/) < p. We observe that at level r"-i(/), either /'-^""'('^(cz) or /'=+'-^""'(')(c2) 
lies in the same component as C2. We consider the two cases: if r"'~^(/) is not a marker, 
then we may proceed two iterates to r"'"'"^(/) keeping the image components distinct. 
If r"~^(Z) is a marker, then r"'~^(/) is marked by J9— (r"~^(/) — r"'(/)); the component 
containing fP~^'^" (')~^"('))(c2) and the component containing C2 at r"~^(/) must have 
distinct preimages at level / which are sent to distinct components of t^{1), one of 
which contains C2, and therefore to distinct components at r"+^(i). We conclude that 
/-r"+i(/) <p. 

Continuing until t°^'^'^^\1) = completes the proof that I < p. □ 

Lemma 2.5. Ifr has period p, and if a level I is marked by k = p — 1, then I < p — 2. 

Proof. Suppose / is marked by p — 1. From Lemma [2.41 / < p — 1. By periodicity, 
r(Z) = 1 — 1. From the admissible r rules, it follows that T{n) = n — 1 for all 1 < n < L 
It follows that n cannot be a marker for any n < l — l. Consequently, level n is marked 
hj I — n for all < n < Z — 1; in particular, / marks level 0. Therefore I 7^ p — 1, 
because p—1 marks level /. 

Proof of Theorem 12.11 Suppose r is periodic with period p. By definition, there 
exists N{t) so that r(n) = n — p for all n > N{t). From Lemma \2.A\ there are no 
marked levels / > p. Therefore, there are no markers at levels / > p + p — 1 = 2p — 1. 
Consequently, r(n + l) = T{n) + 1 for all n > 2p — l, and so we must have r(n) = n—p 
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for all n > 2p—l. If 2p—2 is a marker, then r(2p— 2) = p—1, but this would imply that 
level p — 1 is marked hj p — 1, violating Lemma 12.51 Therefore, r (2p — 2) = p — 2. □ 

Lemma 2.6. Suppose r has length N and t{N) = N — p. Then r extends uniquely 
to a sequence of period p, by setting 

T{n) = n — p 

for all n > N. 

Proof. The existence of the extension follows directly from Theorem I2.2| the unique- 
ness from property (B). □ 

Lemma 2.7. Let r have period p, and suppose T(no) > ''^o — P o-nd T{nQ + 1) = 
uq + 1 — p. Then there exists a marker m < p so that T{m) = uq — p. 

Proof. By periodicity, there is some iterate k so that T^(rio) = UQ—p. By assumption, 
k > 1. Let m = T^~^{no). Because r(no + 1) = uq — p + 1, we have that uq is a 
marker, so m is marked. By Lemma 12.41 then, m < p. We need to show m is also a 
marker. Indeed, r(m + 1) = r(r'^~^(rio) + 1) 7^ ~ P + 1 by property (D). □ 

2.4. Examples/Exceptions. As demonstrated in Theorem 12.11 all periodic tau- 
functions of period p must satisfy r(n) = n — p for all n > 2p — 2. In fact, most 
periodic tau-functions of period p also satisfy r(n) = n — p for all n > 2p — 5. The 
following lemmas provide a complete list of the exceptions. In the lemmas, we express 
the tau-function as a sequence of the form r(l), r(2), r(3), ■ ■ ■ . We remark that these 
lemmas are not used in the algorithm for the Euler characteristic computation, but 
we include them for completeness. 

Lemma 2.8. For each period p > 3, there is a unique periodic tau-function with 
T{n) = n — p for all n > 2p — 2 and T{2p — 3) ^ p — 3. It is given by 

• 0,1,2,-- - ,p-3,0,l,2,--- ,p-2,p-2,p-l,p,--- . 

Proof. By Lemma [2. 7[ there is a marker m < p with r(m) = p — 3. Thus m can only 
be p — 2 or p — 1. Consequently, the tau-function must begin with 0, 1, 2, . . . , (p — 3) 
or with 0, 0, 1, . . . , (p — 3). In the first case. Theorem 12.21 implies that it can only be 
extended as 0, 1, 2, . . . , — 3), with r(2p — 3) = p — 2 and r(2p — 2) = p — 2. In 
the case of 0, 0, 1, 2, . . . , (p — 3), if p is even, then Theorem 12.21 implies the extension 
must be as 0, 0, 1, . . . , (p — 3), 1, 2, . . ., with r(2p — 3) = p — 2, but we cannot extend 
by t{2p — 2) = p — 2. If p is odd, then we must have r(p) = 0, but then T{n) = n —p 
for all n > p. □ 

Lemma 2.9. For each period p > A, the only periodic tau-functions with T{n) = n—p 
for all n > 2p — 3 and T{2p — 4) ^ p — 4 are 

• 0,1,2,--- ,p- 4,0,0,1,2,- ■■ ,p-3,p-3,p-2,p-l,---; 

• 0,1,2, ■■■ ,p-4:,0, 1,2, ■■■ ,p-3,p-3,p-3,p-2,p-l,--- 
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and if p is odd then also 

• 0,0,1,2,-- - ,p-4,l,2,-- - ,p-2,p-3,p-2,p-l,---. 

Proof. By Lemma 12.41 we must have T{2p — 4) = p — 3 or p — 2 or p — 1. Also, by 
Lemma 12.71 level p — 4 is marked by a marker m < p. 

Assume r(2p — 4) = p — 3. Then p — 3 is marked by p — L Periodicity implies 
that p — 4 is marked by 1; that is t{p — 3) = p — 4. The r rules then imply that 
r(n) = n — 1 for all 1 < n < p — 3, so our tau-function begins as 0, 1, 2, . . . , (p — 4). 
Because p — 4 must be marked, Theorem 12.21 implies that t{p — 2) = 0. Theorem 12.21 
then allows for r(p — 1) = or 1. In either case, the tau-function is then uniquely 
determined by Theorem 12.21 and Lemma 12.31 giving the first two possibilities stated 
in the Lemma. 

Now assume r(2j9 — 4) = p — 2. Then level p — 2 is marked by p — 2, so by periodicity, 
we must have t{p — 2) = p — 3 or r(p — 2) = p — A. If r(p — 2) = p — 3, then the 
tau-function begins with 0, 1, 2, ... ,p — 3, but then p — A cannot be marked by a 
marker m < p (contradicting Lemma [2.71) . We must have r(p — 2) = p — 4, and the 
tau-function begins as 0, 0, 1, 2, ... ,p — 4. li p is even, then we can only extend by 
(for p — 4 to be marked), but then r(2p — 4) < p — 3. If p is odd, then we can extend 
by t{p — 1) = 1, and the final tau-function stated in the Lemma is admissible. 

The final possibility is that r(2p— 4) = p—1. The only way to mark p—A by a marker 
m < p is for r(p — 1) = p — 4, so the r sequence begins with 0, T2, rs, 1, 2, . . . ,p — 4, 
for some T2, < 1. Then, as p — 4 is marked by p — 1, we must have r(p) < 2 by 
Theorem 12. 2[ But then r(2p — 4) < p — 1, so its r orbit does not encounter any 
markers larger than 1, and we cannot have T{2p — 3) = p — 3. □ 

Lemma 2.10. For each period p > 5, the only periodic tau-functions withT{n) = n—p 
for all n > 2p — A and r(2p — 5) p — 5 are 

• 0,1,2,--- ,p-5,0,l,2,--- ,p-4:,p-A,p-A,p-3,p-2,---; 

• 0,1,2,--- ,p-5,0,0,l,2,-- - ,p-4,p-4,p-4,p-3,p-2,■■■; 
. 0,1,2,--- ,p-5,0,0,0,l,2,--- ,p-4,p-4,p-3,p-2,■■■; 
. 0,1,2,--- ,p-5,0,l,0,l,2,--- ,p-4,p-4,p-3,p-2,---; 

and if p is odd, then also 

• 0,0,1,2,-- - ,p-5,0,l,2,--- ,p-3,p-4,p-3,p-2,---; 
and if p is even, then also 

• 0,0, 1,2, ■■■ ,p-5, 1,1,2,--- ,p-3,p-4,p-3,p-2, 
and if [p — 1) is divisible by 3, then also 

• 0,1,0,1,2,--- ,p-5,2,3,--- ,p-2,p-4,p-3,p-2,---; 
and if [p — 2) is divisible by 3, then also 

• 0,0,1,1,2,--- ,p-5,2,3,--- ,p-2,p-4,p-3,p-2,---. 



10 



LAURA DE MARCO AND AARON SCHIFF 



Proof. By Lemma \2.7\ p — 5 is marked by a marker < p. By Lemma \2A\ we must 
have t{2p — 5) equal to p — 4, p — 3, p — 2, or p — 1. 

Assume T{2p — 5) = p — 4. Then level p — A is marked by p — 1. Periodicity implies 
that t{p — 4) = p — 3, and therefore that r(n) = n — 1 for all n < p — 4. Therefore, r 
begins with 0, 1, 2, . . . ,p — 5, 0. To reach r(2p — 5) = p — 4, we must have 1 < r(p) < 3, 
allowing only the first four possibilities listed in the Lemma. 

Assume T{2p — 5) = p — 3. We must have t{p — 3) equal to p — 4 or p — 5, by 
periodicity; for p — 5 to be marked, we must have t{p — 3) = p — 5. The tau-function 
begins with 0, 0, 1, 2, . . . ,p — 5. If p is odd, it can be continued by setting r(p — 2) = 0, 
. . . , t{2p — 5) = p — 3, and r(2p — 4) = p — 4. If p is even, then t{p — 2) = 1, so we 
can take r(p — 1) = 1 to allow for r(2p — 5) = p — 3. 

A similar argument handles the case of T{2p — 5) = p — 2. □ 



3. Generating periodic tau-functions 

The goal is to generate a list of all tau-functions of period p. For the later steps in 
the algorithm, we need the data of the tau-functions themselves, not only the total 
number. 

There is a unique tau-function of period p = 1, given by r(n) = n — 1 for all n > 1. 
For small periods, say period p < 10, there are few periodic tau-functions. Applying 
Theorem 12. 2[ we can generate all tau-functions to length 2p — 2. From Theorem 12. ![ 
the equality r(2p — 2) = p — 2 holds if and only if this tau-function extends to a 
sequence of period p; further, the extension is uniquely determined. For example, 
the total number of tau-functions of length 8 (= 2p — 2 for p = 5) is only 144, so 
the computation time and memory usage are negligible for period p = 5 [DS]. As 
the period grows, the total number of tau-functions of length 2p — 2 grows fast; it is 
probably larger than 4^"^. We use the Lemmas of the previous section to reduce our 
computational requirements. 

3.1. Algorithm. The algorithm proceeds as follows. Fix p > 1. We generate a list 
called Periodic containing all tau-functions of period p. For the induction step, we 
generate a list called Continue. 

Initialization. Generate all tau-functions to length n = p, following 
Theorem 12.21 If r(p) = 0, include in Periodic. If r has no markers, 
then discard. Otherwise, include in Continue. 

Extension to length n + 1 and test for periodicity. Choose r 
from the list Continue. Let n be its length; by construction, r(n) > 
n—p. Determine values /q, /i, . . . (as appearing in Theorem 12. 2 [ setting 
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N = n) subject to the extra condition n — p < li < Iq = riji). For 
each such /j, we consider the admissible extension of r, defined by 

r(n + 1) = li + l. 

li li = n — p, then include the extended r in Periodic; by Lemma 12.61 
this T uniquely determines a periodic tau-function. 
If n < 2p — 3 and if k > n — p and if 

max{r(m) : m < p, T{m + 1) < r(m)} > n — p, 

then include in Continue; this r is a candidate to have a periodic 
extension, as it satisfies the necessary conditions of Lemmas 12.31 and 
12.71 and Theorem 12. 1[ Otherwise, discard. Repeat the induction step 
until Continue is empty. 

3.2. Details. In Tables [2] and [3l we include the particulars of our computation for 
generating all periodic tau-functions of periods 10 and 20. Following the algorithm 
above, we show the number of tau-functions in the lists Periodic and Continue as we 
increase the length of the tau-functions. 



Length 


Periodic 


Discard 


Continue 


10 


205 


1 


435 


11 


201 


242 


506 


12 


139 


567 


479 


13 


57 


780 


279 


14 


26 


497 


134 


15 


12 


251 


61 


16 


6 


122 


21 


17 


2 


43 


6 


18 


1 


13 






Table 2. Period 10 details: generating the 649 tau-functions of period 
10 from a total of 279,415 tau-functions of length 18. Final data file 
size = 7.8 KB, peak disk usage = 18 KB. 

4. Topological conjugacy classes of basins 

In this section, we describe the algorithm to compute the number Top(r) of topo- 
logical conjugacy classes of basins (/, X(/)) with a given tau-function r. It is proved 
in |DP] that Top(r) can be computed as 

Top(r) = Spines(r) ■TF(r), 
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Length 


Periodic 


Discard 


Continue 


20 


449308 


1 


848362 


21 


319756 


528624 


1055320 


22 


389254 


1059653 


1116657 


23 


114128 


1523035 


978211 


24 


41925 


1646071 


674730 


25 


17081 


1299907 


391444 


26 


8896 


800601 


196937 


27 


4138 


403194 


93346 


28 


1898 


192799 


44601 


29 


978 


92478 


20839 


30 


475 


43078 


9636 


31 


217 


20028 


4571 


32 


113 


9623 


2054 


33 


52 


4309 


932 


34 


24 


2004 


414 


35 


12 


901 


169 


36 


6 


373 


57 


37 


2 


137 


18 


38 


1 


41 






Table 3. Period 20 details: generating the 1,348,264 tau- functions of 
period p = 20 from a total of about 1.5 trillion tau-functions of length 
2p - 2 = 38. Final data file size = 29 MB, peak disk usage = 74 MB. 

where Spines(r) is the number of pictographs (or truncated spines) associated to r and 
TF(r) is the associated twist factor. We include here the steps to compute Spines(r) 
and TF(r). These details already appeared in |DS] . 

The twist factor TF(r) is denoted by Top(P) in [DP], the number of conjugacy 
classes of basins with pictograph V, for any pictograph with tau-function r. Indeed, 
it is easy to see that any pictograph with a period tau-function will have only finitely 
many marked levels, thus satisfying the hypotheses of |DPt Theorem 9.1]; further, 
it is stated there that the computation of Top{T>) depends only on the underlying 
tau-function. 

4.1. Computing the number of pictographs. Fix an admissible tau-function r 
of length N. As in §2.H the markers of r are the integers 

{m e {1, . . . ,N - 1} : T{m + 1) < r(m) + 1}. 

The marked levels of r are all integers in the forward orbits of the markers: 

{/ > : / = T'^{m) for marker m and n > 0} U {0}; 
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we say is marked even if there are no markers. 

As in Theorem 12.21 we let k be the number of markers which appear in the orbit 

^ t{N) ^ ...^ T^'^^WiN) = 0. 

Label these k markers hj l[, I2, . . . , I'l^ so that 

AT = > /; > > . . . > 4 > 0. 

For each < i < k, let U = r(/-) so that 

t{N) = lo>h> ■■■>lk>0. 

For each < i < k, define rii by the condition that 

and define so that T"'*(/fc) = 0. (The Ui are called special orders in the program.) 
For < i < j < + 1, we set 

6H ,) = l^ ifr(/: + l) = /, + l 
\ otherwise 

where by convention we take /fc+i = —1- Note that r(Z^ + 1) = for every r, so 
S{k,k + 1) = 1. 

The symmetry of r is 

s = min{n > : r"(/o) is a marked level}. 

Note that s < uq. To each admissible choice for t{N + 1) (from Theorem 12.21) we 
define the (A^ + l)-th spine factor of r. If t{N + 1) = /j + 1 with i > 0, we set 

SF(r, N+l):= 2"«-^(2"i (2"^(- ■ ■ (2"-i -5{i-l,i)) ) - S{2, i)) - 5(1, i)); 

as above, we take Ik+i = — 1. If t{N + 1) = Iq + 1 = t{N) + 1, we set 

SF(r,A^ + 1) = 1 

The number of pictographs (or equivalently, truncated spines) associated to a tau- 
function is computed inductively on the length. 

Proposition 4.1. Let t be a periodic tau-function of period p. The number of pic- 
tographs with tau-function r is given by 

N 

Spines(r) = nSF(r,j). 
for any choice of N with t{N) = N — p. 
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Proof. That Spines(r) is the product of spine factors is deduced in |DSj . It remains 
to show that the computation terminates at a finite N when r is periodic. From the 
definition of the spine factor, it is equal to 1 whenever t{N + 1) = t{N) + 1. For 
periodic taus, this will be the case for all N sufficiently large. Recall from Lemma [273] 
that once we find one N with t{N) = N — p, this equality will hold for all n > N. □ 

4.2. Computing the twist factor. Fix an admissible tau-function r of length G 
NU {oo} with finitely many marked levels. For each n < N, the order of n was defined 
in §2.2( it satisfies r™'^*^")(n) = 0. For each marked level / > 0, compute 

I 

mod(/) = ^2-°'^'i« 

i=l 

and 

t{l) = min{n > : nmod(/) G N}. 
We define the twist period T{t) by 

(4.1) T{t) = max{)f:(/) : / is a marked level} 

or set T{t) = 1 if r has no non-zero marked levels. 

Let L{t) be the number of non-zero marked levels. The twist factor is defined by 

2L(t) 

^ nry 

Theorem 9.1 of |DP] states that the number of topological conjugacy classes of basins 
associated to a given pictograph with tau-function r is equal to TF(r). 

4.3. The significance of the twist factor. We include a few words here to explain 
the meaning of the values appearing in §4.21 to define the twist factor. These play a 
role in the explanations of ^ 

The quasiconformal deformations of the basin of infinity of a polynomial / have a 
natural decomposition into twisting and stretching factors; see |McS] or the summary 
in |DP] . Let / = f^a,v) be a cubic polynomial with periodic tau-function r. Let be 
its escape-rate function. Recall that —a is the critical point that escapes to infinity. 
The fundamental annulus of / is the domain 

A{f) = {ze€: Gf{-a) < Gf{z) < 3Gf{-a)}. 

Viewing the basin of infinity X{f) as an abstract Riemann surface, a full Dehn twist 
in A{f) induces the hemidromy action described in |BH] : see also [Br] for an accessible 
summary. 

The twist period T(r) is the least power of a full Dehn twist in the fundamental 
annulus that lies in the mapping class group of /. To compute T(r), we determine 
the induced amount of twisting in any image or preimage of A{f) under the action 
of /. Following the descriptions in [Br] and |BH] . it suffices to compute the relative 
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moduli of these annuli lying between the two critical points; the relative modulus of 
an annulus A is the ratio mod(A)/ mod{A{f)). The value mod(/) computes exactly 
these sums of relative moduli down to the l-th marked level. 

The twist factor is the ingredient emphasized in |DP] . By measuring twist periods 
against the total number of ways to produce basins (/, X{f)) from a given pictograph, 
the discrepancy amounts to the twist factor TF(r). 



5. Escape regions 

In this section we explain the final steps of the algorithm, incorporating the com- 
putations described in the previous section. 



5.1. The moduli space. As discussed in [Mi] , there is a natural involution on the 
space V^"^, given by 

/(a, v) = (—a, —v) 

induced by the conjugation of f(a,v) by z i-t- —z. Thus there is a degree 2 projection 

to the moduli space of critically-marked cubic polynomials. The action of / preserves 
the curve Sp, defining a curve Sp/ 1 C A^g™. 

5.2. Escape regions and multiplicity. As introduced in |Mi] and |BKM] . an escape 
region of Sp is a connected component of 

{f{a,v) e Sp : ^ oo as n ^ oo}. 

That is, it consists of maps with one periodic critical point (at +a) and one escaping 
critical point (at —a). 

It follows from the general theory on stability that all polynomials in a given escape 
region E are topologically conjugate on C, as described in |McS] . In this special 
setting, though, it can be seen directly from a canonical parameterization of E. It is 
shown in |Mit Lemma 5.6] that each escape region E is conformally a punctured disk, 
canonically identified with an m-fold cover of a punctured disk, for some positive 
integer m = m{E). This number m{E) is called the multiplicity of E. 

The covering map of degree m{E) is defined by the assignment 

where (p(a,v) defines the uniformizing Bottcher coordinates near infinity for f(a,v): 
where f(a,v){f(a,v)iz)) = {ip(^a,v){z)Y , unique if chosen to satisfy v^'(c)o) = 1. The 
point 2a is the cocritical point for —a, so /(a^„)(2a) = f(a,v){—ci)- In particular, the 
twisting deformation on the basin of infinity induces the change in angular coordinate 
on E. In fact, the external angle of 2a is increased by vr under a full Dehn twist in 
the fundamental annulus of f(a,v)', thus, 2m{E) full twists closes a loop in E. 
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Lemma 5.1. Fix an escape region E and let r he the tau-function of any f G E. 
The multiplicity is given by 



Proof. Each escape region E projects to an escape region E/I m the curve Sp/I C 
JVi"^. By definition of the twist period, T(r) full twists in a fundamental annulus are 
required to induce a closed loop in E/I. But E/I is doubly covered by a single escape 
region E if and only if /(a,i,) and f(-a-v) are equivalent under a twist deformation, if 
and only if we have T(r) = 1. In this case of T(r) = 1, two full twists are required 
to close a loop in E, corresponding to an argument increase of 27r for the cocritical 
point 2a. Therefore m{E) = 1. On the other hand, if T(r) > 1, then each escape 
region E projects bijectively to E/I; thus 2m{E) = T{t). □ 

5.3. Hybrid classes. For any polynomial / in an escape region E in Sp, the associ- 
ated tau-function will have period k for some k dividing p. A restriction of the iterate 
to a certain neighborhood of +a will then define a quadratic polynomial-like map. 
We refer to |DH] for background information. In this context, it is important to know 
that the conformal conjugacy class of / is uniquely determined by the conformal con- 
jugacy class of its basin of infinity (/, X{f)) and the hybrid class of its polynomial-like 
restriction |BH] . See also |BKMt Theorem 3.9, Corollary 3.10]. 
We will use the following consequence of the general theory: 

Proposition 5.2. An escape region E/I in Sp/I is uniquely determined by 

(1) an integer k dividing p with 1 < k < p; 

(2) a topological conjugacy class of basin dynamics {f,X{f)) with a critical end 
of period k; and 

(3) a point in the Mandelbrot set corresponding to a center of period exactly p/k. 

A center of period n in the Mandelbrot set is a solution c to the equation /"(O) = 
where fc{z) = + c. The center c has period exactly n if n is the smallest positive 
integer for which the equality /"(O) = holds. The number V2{n) of centers of period 
exactly n is easily computable by the following relation: 



Combining the above results, we deduce the following: 

Theorem 5.3. For any tau-function r with period k dividing p, the number of escape 
regions in Sp with tau-function r is 




where T{t) is the twist period computed in §^.4 



2 



in— 1 



Ml) 



q\n, l<9<n 



Ends{T,p) = 



U2{p/k) Spines(r) TF(r) «/T(r) = 1 
2z/2(p/A;) Spines(r) TF(r) «/r(r) > 1 
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where T{t) is the twist period, Spines(r) is the number of pictographs, and TF(r) is 
the twist factor ofr. The total number of escape regions in Sp is therefore 

Np = J2 Yl Ends{T,p) 

k\p per{T)=k 

In particular, in the case of k = p, Ends{T,p) is the number of "central ends" of 
r, coinciding with the number of all escape regions of Sp with tau-function r and 
hybrid class z"^. The sum of Ends{T,p) over all taus with period p is shown in Table 
[H The sum of Ends{T,p) over all taus with period dividing p is the total number Np 
of escape regions in Sp. 

Proof. Fix r of period k dividing p. From the arguments of §11 there are Spines (r) TF(r) 
topological conjugacy classes of basins {f,X{f)) of cubic polynomials with tau- 
function r. Applying Proposition l5.2l there are consequently voiplk) Spines(r) TF(r) 
escape regions E/I in Sp/I C Mf^. If T(r) = 1, then exactly as in the proof 
of Lemma 15.11 there is a unique escape region E in Sp mapped to each E/I. If 
T(r) > 1, there are exactly two escape regions mapped to each E/I. □ 

5.4. Testing the computation. We conclude with an explanation of the test of our 
computation against the degree of Sp. 

The multiplicity of an escape region E in Sp coincides with the number of intersec- 
tion points of E with any line in Vf^ of the form {a = a^] for any Oq of sufficiently 
large modulus. Therefore, the degree dp of the curve Sp must satisfy 

dp = ^m{E), 

E 

summing over all escape regions E of Sp. The degree dp is easily computed, as it 
satisfies: 

where the sum is taken over all q dividing p with 1 < q <p. As established by Lemma 
15. H the value m{E) depends only on the tau-function for the escape region E, so we 
may define 

m(r) := m{E) 

for any escape region E associated to tau-function r. 

Our algorithm determines the value Ends{T,p) for every tau-function of period k 
dividing p; the ingredients are listed in Theorem 15. 3[ We can therefore check our 
computation by assuring equahty of 

m(r) Ends{T,p) = dp, 

T 

summing over all tau-functions r of periods dividing p. 
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